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This manuscript was not published eventually, for two 
reasons: 
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• The findings were too similar to the localization of 
wrinkles obtained earlier for a sheet on an elastic foun- 
dation. See G. W. Hunt, M. K. Wadee, and N. Shia- 
colas, J. Appl. Mech. 60, 1033 (1993). Note, however, 
that the two problems are not equivalent. For an elas- 
tic foundation the second term in the energy density, 
eq. (1), is (K/2)h 2 rather than (K/2)h 2 cos 0. 

• We have found the exact solution to the problem pre- 
sented here. See arXiv:1107.5505. 

In addition, B. Audoly has followed this preprint to 
^ ' show, using an amplitude-equation approach, that the vari- 
^ , ational Ansatz used here is indeed the profile minimizing 
r ^ ■ the energy to quartic terms in the height amplitude. See 
B. Audoly, Phys. Rev. E 84, 011605 (2011). This article 
also presents a unified pcrturbative treatment for several 
types of substrate, including the fluid substrate and the 
linear elastic foundation, 
i— —i ■ If you find the manuscript useful nonetheless, and wish 

' to cite it, please use larXiv:1009. 24871 
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Abstract. We analyze the buckling of a rigid thin membrane floating on a dense fluid substrate. The 
interplay of curvature and substrate energy is known to create wrinkling at a characteristic wavelength A, 
which localizes into a fold at sufficient buckling displacement A. By analyzing the regime 4 C A, we show 
that wrinkles are unstable to localized folding for arbitrarily small A. After observing that evanescent waves 
at the boundaries can be energetically favored over uniform wrinkles, we construct a localized Ansatz state 
far from boundaries that is also energetically favored. The resulting surface pressure P in conventional 
units is 2 — (it 2 /4)(A/ A) 2 , in entire agreement with previous numerical results. The decay length of the 
amplitude is k~ - 1 = (2/tt 2 )A 2 /A This case illustrates how a leading-order energy expression suggested by 
the infinitesimal displacement can give a qualitatively wrong configuration. 



PACS. XX.XX.XX No PACS code given 



1 Introduction 

Recent work has revealed a wealth of spontaneous spatial 
structures in thin elastic sheets, owing to their easy de- 
formability by bending [^^Pl[51|o1[7l^[Tu1P1P^[T^[T^ 
[Po] . Special interest has focusses on the recently-proposed 
wrinkle-to-fold transition in compressed thin sheets float- 
ing on a dense liquid [TO]. Such composite structures oc- 
cur widely in industrial coatings and biological tissues. 
The sheet can be a 50-nm-thick metal film j8j or a two- 
nanomcter-thick lipid monolayer [SJ. 

When compressed to the point of buckling, the incipi- 
ent deformation is an undulation or wrinkle whose wave- 
length A is a combination of the bending stiffness and 
the substrate rigidity, as specified below. A compressional 
displacement A accompanies this buckling. Buckling oc- 
curs when the compressional force per unit length P ex- 
ceeds a threshold. In the conventional units defined below, 
this threshold P is equal to 2. As with the classic Euler 
buckling of a compressed rod |16j . the buckled structure 
is unstable under this force: larger displacement creates 
more bending deflection, which requires a smaller pres- 
sure. Thus P decreases as A increases. 

For sufficiently large A of order A, the weak-deflection 
expression for the energy becomes inadequate, and the se- 
lected structure is no longer a uniform undulation but 
a localized fold. Reference [TUJ explored the transition 
from the wrinkled to the folded state numerically and 
experimentally, concluding that the transition occurs for 
A ~ A/ 3. They determined numerically that the dimen- 
sionless force P is well represented as 2 — 2.47(Z\/A) 2 . As A 



decreases, the buckled region was observed to be less and 
less localized, and to span more and more wavelengths, 
eventually spanning the system. 

These observations led us to consider the regime of 
indefinitely small displacements A <C A, in order to dis- 
cern how the corrections to the small-deformation energy 
emerge. As explained below, we found that these correc- 
tions remain important even for arbitrarily small displace- 
ments. The resulting buckling pattern is always localized, 
though its spatial extent diverges as the displacement A 
decreases to zero. Since the regime of small displacements 
is analytically tractable, we could obtain a closed-form ex- 
pression for the P(A) function: P = 2 - (vr 2 /4)(Z\/A) 2 + 
C((Z\/A) 4 ). This agrees completely with the numerical 
value cited above. We also obtain a closed-form expres- 
sion for the buckling profile in this regime. Though this 
form is obtained by a variational Ansatz, we argue below 
that it approaches the exact profile in the limit of small 
displacement. 



2 Model 

We consider a two-dimensional problem of a thin incom- 
pressible elastic sheet of length L, width W, and bending 
modulus B. The sheet is assumed to deform in the xz 
plane while remaining uniform along the y axis. Being 
incompressible, its total length L is fixed, and its configu- 
ration, including the total displacement along the x axis, 
A, is fully accounted for by a height profile, h(s), as a 
function of arclength s E [-L/2, L/2]. The region z < h 
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is occupied by a liquid of mass density p — K/g (g being 
the gravitational acceleration). See fig.ffl 




Fig. 1. Schematic view of the system and its parametrization. 

For a certain configuration, h(s), the energy of the 
sheet E and the displacement A are given by 

L/2 



E = W I ds ( + ]-Kh 2 cosi 

J-L/2 V 2 2 
r-L/2 

A 

I -L/2 



,L/2 

/ ds (1 — cos ( 

J-L/2 



(1) 



where <j>(s) is the angle between the local tangent to the 
sheet and the x axis at arclength s (cf. fig. [T]), and a dot 
denotes a derivative with respect to s. It is helpful to work 
in units where B = K = 1, i.e., rescale the energy by B 
and all lengths by (B/K) 1 ^ 4 . The pressure P is then scaled 
by {BKfl 2 . We shall also let W = 1 in what follows. 
Using these definitions and the geometrical relation h = 
sin <fi, we rewrite eq. ([!]) as 



E 



A = 



L/2 



ds 



-L/2 
L/2 



L/2 



1-h 2 



1-h 2 



1 - \ 1 - h 2 . 



( 2 ) 



Our aim is to find the height profile, h(s), which mini- 
mizes the energy of the sheet for a given displacement A, 
E = E*(A); from this we may find the resulting pressure- 
displacement relation, P(A) — d_E*/dA 



3 Extended wrinkling 

Let us substitute in eq. ([2]) a uniformly wrinkled height 
profile of amplitude A and wavenumber k, 



h(s) = A cos(fcs). 



(3) 



The extensivity of the energy implies that E = L£(5) with 
8 — A/L. To second order in A 2 we get 

£ = \{k 4 + 1)A 2 + ^k 2 {2k 4 -1)A 4 , (4) 



S = Ik 2 A 2 + -k 4 A 4 . 
4 64 



(5) 



We first examine the leading order. Substituting eq. 
© in eq. g]) yields £{S) = (k 2 + k~ 2 )5. For any 6 > 
the energy has a minimum of £* — 25 for k = 1. This 
marks the linear wrinkling instability [17,3,18 , with a fi- 
nite wavenumber k c ~ 1 [wavelength A = 2tt(B / K) 1 / 4 ), 
squared amplitude A 2 = AS, and energy E* = 2A. The 
resulting critical pressure is P c = dE* /dA = 2. 

Repeating the same procedure in the next order, we 
find from eqs. dl]) and (0 



£ (5) = (k 2 + k- 2 )8 + (k 2 - 5AT 2 )<5 2 /4. 



(0) 



This expression is still minimized by k = k c = 1, yielding 
£*{5) =2<5-(5 2 , i.e., 



E*(A) = 2A-A 2 /L 
P[A) = 2- 2A/L. 



(7) 



4 Localized folding 



While the solution to the leading-order problem must be of 
the wrinkling form [eq. ([5])]. such an extended profile does 
not necessarily minimize the anharmonic energy. Indeed, 
the next-leading energy correction in eq. (T^Tj) vanishes as 
L — > oo ; the wrinkle profile gains no energetic benefit from 
this correction in the infinite system. 

To see that a localized solution should be favorable, 
let us fix the pressure (rather than the displacement) at a 
value slightly below threshold, P = 2 — p with < p <£L 1, 
and re-examine the leading-order problem. The energy to 
be minimized in this case is Q = £ — P5 ~ A 2 (k 4 — Pk 2 + 
l)/4, where we have used the leading order in eqs. (HJ) 
and ([5]). The most favorable wrinkling profile has k 2 = 1 
and Q = (A 2 /A)p > 0, which is larger than Q = of the 
flat state. Thus, the fiat state is linearly stable against 
wrinkling for P < P c . Nonetheless, we can get a lower 
energy if we allow for a complex wavenumber. We may 
obtain Q = by taking 



k ~ 1 ± in, k = y/p/2. 



(8) 



Hence, according to the leading-order theory, if there were 
a boundary at s = 0, an 'evanescent wave' adjacent to the 
boundary of the form h(s > 0) = Ae ls ~ KS , could be stabi- 
lized for P < P c . For such a fixed pressure the flat state is 
linearly stable against wrinkling, and a 'propagating wave' 
can be stabilized strictly at P = P c . 

We are interested, however, in profiles which are lo- 
calized far away from boundaries. Accordingly, we seek a 
wave-packet profile with a spectral content similar to the 
evanescent-wave profile above. One such profile is 



h(s) = A 



cos(ks) 
cosh(Ks) ' 



(9) 



where the decay coefficient n is used as a variational pa- 
rameter, and we take k — k c = 1. (Letting k be another 
variational parameter does not change the results.) We 
assume the limit L _1 « k C I; = 1, where the wavy pro- 
file is 'underdamped' but decays sufficiently fast for the 
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system size to be taken as infinite. The extensivity of the 
energy again implies E = k~ 1 £(5), where now S = kA. As 
in the preceding section we perform the calculation to or- 
der S 2 ; yet, unlike the extended- wrinkling case, the result 
will be of higher order in A, since the domain size itself, 
ft -1 , depends on A. 

Substituting eq. ([9]) in eq. ([2]), we obtain within the 
assumed order of approximation, 

£ = (l + fi 2 )A 2 + ii 4 (10) 

6 = 1(3 + k 2 )A 2 + ±;A\ (11) 
6 16 

Solving eq. ([II]) for A 2 , 

A 2 = 2(1- k 2 /3)5-5 2 /2, (12) 

and substituting it back in eq. (|10l) . we get £ (S) = 2(1 + 
2k 2 /3)<5 — S 2 /3. To properly minimize the energy with re- 
spect to « for fixed A, we rewrite this result as 

E = 2(1 + 2k 2 /3)A - nA 2 /3. 

As expected, to the leading order the energy is mini- 
mum for k = and coincides with that of the wrinkling 
case, E* = 2A. However, with the anharmonic term the 
energy is minimized by k = A/ 8, yielding 

A 3 

E*(A) = 2A- — 

P(A) = 2-^r. (13) 

Note that the result k = A/8 = \J2 — P/2 is identical to 
the one obtained for the evanescent wave from the leading- 
order problem, eq. ©. 

Finally, substituting k = A/8 back in eqs. (|12p and 
([5]), we obtain for the localized profile, 

, , A( 7A 2 \ coss 

H(S)= 2 I 1 ' 384 J cosh[(Z\/8) S ]- (14) 

Equations (|13|) and (fT4)l are our central results. Com- 
paring the energies of the extended and localized profiles, 
eqs. ([7]) and (fTB"]) . we find that the localized fold is sta- 
bilized for A > A c = 48/L. Thus, the wrinkle-to-fold 
transition presented in ref. |10) is a finite-size effect; in 
the limit L — > oo the deformation of the sheet is always 
localized in a folded domain of negligible width compared 
to the system size, kT x = 8/ A [or, in dimensional terms, 
kT 1 = (2 / tt 2 )\ 2 / A]. The parabolic pressure-displacement 
relation of eq. (|13p is rewritten in dimensional terms as 

p -o «y*y 

(BK) 1 / 2 4 V a; ' 

which agrees nicely with the numerical result of ref. [101 , 
P/(BKf/ 2 ~2-2A7(A/\) 2 . 

From eqs. and (fTB"]) one might conclude that the 
appearance of the localized state at A = A c is a first-order 



transition with discontinuous jumps in k and P. This is 
an artifact, arising from the assumption k" 1 -C L, which 
breaks down for A <~ k ~ As shown in the Appendix, 
a more careful analysis for such small displacements yields 
a second-order transition, with k = for Z\ < Z\ c and 
n ~ (A — Ac) 1 / 2 for zi > A c , leading to a discontinuous 
jump in dP/dZ\. 



5 Numerical solution 

The calculations given in sects. [3] and 2] are valid suffi- 
ciently close to the instability. To examine the deformation 
of the sheet at larger displacements we derive the Euler- 
Lagrange equation, corresponding to the minimization of 
E for a given A, and solve it numerically. 

We would like to find the profile h(s) that minimizes 
the elastic energy E under the constraint of fixed displace- 
ment A, where E[h(s)] and Z\[/i(s)] have been defined in 
eq. ([2]). The following equation is obtained from calculat- 
ing the variation of G = E — PA with respect to h(s) and 
setting it to zero: 

(1 - h 2 ) 2 "h + 4(1 - h 2 )hhh + (1 + 3h 2 )h 3 

+ (h 2 /2 + P)(l - h 2 f/ 2 h + h(l - h 2 f' 2 = 0. (15) 

This equation needs to be supplemented by four boundary 
conditions, for example, the hinge conditions h(±L/2) = 
and h(±L/2) = 0. 

In the following examples we look for even profiles, 
h(—s) = h(s), over a system length of L = 2\-n (i.e., 
containing 10 ^ w rinkling wavelengths). For a given P < 
P c = 2, eq. (|15[) is solved numerically over the interval 
s € (— L/2,0) using the following boundary conditions: 
h(-L/2) = 0, h(-L/2) = 0, h(Q) = 0, h(0) = h . We 
shoot for the value of ho so as to make | ft,' (0)1 vanish 
within the required accuracy. Once h(s) is found, the dis- 
placement A(P) is calculated using eq. ©. The results 
are presented in fig. [51 where they are compared with the 
analytical ones, eqs. (fT3)l and ([T4)l . 

The discrepancy between the numerically determined 
height profile and the theoretical one is of the order of a 
few percent, up to about 10%. It increases with A as ex- 
pected [fig. (UJa)]. The predicted parabolic law for P(A), 
however, agrees with the numerics to a much greater ex- 
tent [fig. [2Kb)]. The agreement is remarkable for large dis- 
placements, where the perturbative theory is not expected 
to hold at all. Since we have carried out the analytical cal- 
culation to order S 2 , we should have expected a correction 
to E of order k~ 1 5 3 ~ A 5 , leading to a correction of order 
A 4 in P. Even for nearly vertical deflections (right side of 
fig. [3]) the discrepancy is consistent with numerical error. 
The large discrepancies at small A come from the finite-L 
effect discussed above. This effect should become exponen- 
tially small in kL ~ LA, which explains the sharp decay 
of the discrepancy a a with A and its strong dependence 
on L. [See fig. EJb) inset.] We have replotted in fig. [3] the 
data of that inset as a function of kL = (A/8)L on a semi- 
logarithmic scale, confirming the exponential decay of the 
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Fig. 2. (a) Height profiles obtained from the numerical solution (black solid curves) and the approximate theory [eq. (|14p . 
blue dashed curves]. Different curves correspond (from top to bottom) to decreasing pressure (as indicated) and increasing 
displacement. The length is L — 2ln. Consecutive curves are vertically shifted by 1.5 for clarity. Note that the profiles are 
presented in terms of the material coordinate s rather than the projected coordinate x. (b) Pressure as a function of displacement. 
The open circles are the numerical values obtained from the profiles shown in (a). The solid curve is the theoretical prediction 
[eq. (1131) ]. The inset shows the normalized discrepancy a a between our numerical A(P) and that of eq. (| 13p : specifically, 
a a = (4mm — Ath)/A num . Solid circles are for L — 21n; open squares are for L = ll7r. 



discrepancy with kL. We note that the numerical results 
presented in ref. [TU] seem to indicate that the parabolic 
pressure-displacement relation remains accurate up to the 
point of self-contact. All of the above suggests that P(A) 
of eq. (|13p might, in fact, be the exact relation in the limit 
kL > 1. 



6 Discussion 

In this section we argue that for the limit of small am- 
plitude our Ansatz should become exact. We discuss the 
meaning of the simple large-amplitude behavior noted above 
We comment on the relevance of our result for recent 
Langmuir buckling observations. We draw parallels to other 
continuous instabilities that are not predicted by lowest- 
order perturbation of the energy function. 

Figure [2] shows a near match between the numerically 
determined profile and the Ansatz of eq. (0) ■ The agree- 
ment becomes better as the threshold is approached. In 
this limit our asymptotic 0(A 4 ) energy becomes accurate. 
Further, the spectral profile of eq. © approaches that of 
the evanescent wave that minimizes this asymptotic en- 
ergy. This profile has a peak at k = 1 of width k. The 




KL = LA/8 



Fig. 3. Discrepancy between theory and numerics as a function 
of nL = (A/8)L for L — 21n (solid circles) and L = Htt (open 
squares) on a semi- logarithmic scale. The fitted line has a slope 
of -1.03 ± 0.04. 

profile of eq. ([§]) has a spectral weight that falls off ex- 
ponentially for k — I k. We expect that any Ansatz 
with such a spectral distribution would approach the ac- 
tual profile near the threshold. However, not any profile 
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with the proper width and peak position is adequate. One 
such profile is h(s) — A sins e _K ' s '. This profile gives the 
same qualitative behavior of the energy and pressure for 
a given displacement A as eq. ©, yet, here the spectral 
weight falls off as a power of k — 1, rather than the ex- 
ponential falloff of eq. ©• Accordingly, we found that its 
energy is fractionally larger than that of eq. (|9]). 

Our Ansatz is not unique. For example, the alternative 
form 

h(s) = Acos(ks)/ cosh 2 (ns/2) 

has the same spectral features noted for eq. (J9]). However, 
neither profile is exact; thus when expanded in powers of 
A, both must differ from the true behavior in some order. 
For example, the (7/384)Z\ 2 correction in eq. (fT4]) depends 
on our choice of Ansatz. 

The theory predicts a quadratic initial dependence of 
pressure on displacement, as reported already in ref. |10j . 
As noted above, the quadratic dependence remains accu- 
rate far beyond the small-displacement regime where it 
was derived. This accuracy suggests that the system may 
be described in a way that is manifestly homogeneous in 
A. However, we have not been able to find a simple theory 
of this type. 

The main experimental implication of the theory above 
is that one should not expect to see a transition from a 
wrinkled to a folded state in a large experimental system. 
For overall displacements A of order A this finding is im- 
plicit in the work of ref. |10j . The strain and the wrinkle 
amplitude at the threshold of their fold behavior goes to 
zero as the system size L goes to infinity. Nevertheless, 
our work shows that the transition to a strongly folded 
state, presented in ref. |10) . is a smooth crossover, not a 
true transition. Even when the overall displacement A is 
indefinitely smaller than the wrinkling wavelength A, the 
uniformly wrinkled state is unstable against a state of lo- 
calized deformation. 

Furthermore, if the assumption of incompressibility 
employed throughout our analysis is relaxed, another non- 
zero critical displacement appears, where a transition from 
the flat to a wrinkled state takes place. This buckling oc- 
curs at a displacement A^ where the in-plane compres- 
sional stress Y(A\ 3 /L) (Y being the two-dimensional com- 
pression modulus of the sheet), exceeds the buckling stress 
of the sheet 2(BK) 1 ^ 2 . Now, if the system size L is suf- 
ficiently large, we will have A c < A\>, and the regime of 
uniform wrinkling will disappear altogether. The condi- 
tion for such a direct flat-to-fold transition [TB] is L > 
(y/iQV2 ^ (Y/B)^ 2 X 2 ~ X 2 /t, where we have used the 
scaling of the compression and bending moduli with the 
sheet thickness t: Y ~ t and B ~ t 3 . The inclusion of 
compressibility renders the transition first-order [16] . in 
contrast to the continuous incompressible case treated by 
ref. [TU] and in the present work. 

The wrinkled state has not been observed in macro- 
scopic Langmuir monolayers [S]. We believe that the the- 
ory above explains this lack of observation. In a typical 
Langmuir- monolayer experiment L ~ 10 cm and A ~ 1 
fan. Hence, a wrinkled state should get localized already 
for extremely small displacements, A > A c ~ \ 2 /L ~ 



0.1 A. The thickness of the monolayer is of 1 nm scale, im- 
plying the lack of extended wrinkling for L > 1 cm. Thus, 
we expect that in such experiments there is a direct transi- 
tion from flat to folded configurations. Because of the rele- 
vance of compressibility in this experiment, we expect the 
folds to form via discontinuous, abrupt events, as observed 
[5]. The wrinkled state has been reported experimentally 
10, 13] in systems that are not extremely large compared 
to A. As discussed in the Appendix, the wrinkled state can 
be stable in systems where L < A 2 /A. Other experiments 
exhibiting wrinkling [7] have strongly constrained lateral 
boundaries that likely form a barrier against folding. 

The wrinkle-to-fold transition has been applied to de- 
scribe not only fluid substrates but also elastic substrates 
of stiffness K. In the lowest-order approximation, these 
two situations are indeed equivalent. However, in the next- 
leading order treated here, the elastic system becomes dif- 
ferent from the fluid system because an elastic substrate 
then exerts shear stresses on the sheet, while the fluid 
substrate does not. It may well be that the wrinkled state 
is more stable on an elastic substrate than on the liquid 
substrate treated above. 

It is interesting to compare the present instability with 
that seen in an ordinary first-order transition near the 
spinodal line [19] , In this case as in ours, the leading- 
order contribution to the energy is stabilizing — favor- 
ing zero amplitude of the new phase. In both cases the 
next leading order is destabilizing. In both cases one may 
choose an intensive variable (i.e., pressure) as the control 
parameter, or an extensive one (i.e., volume or our dis- 
placement A). In the phase transition case, the current 
phase is metastable. A nonzero threshold volume of ma- 
terial must nucleate the new phase in order for the tran- 
sition to occur. In the folding transition, there is likewise 
a critical displacement required for any pressure P be- 
low the threshold. As the control parameter moves to the 
spinodal line, the new phase dictated by the lowest-order 
energy begins to grow spontaneously. The counterpart of 
this behavior in our system would be for a wrinkled phase 
to grow in amplitude. However, the above analysis indi- 
cates a different behavior for the case of folding. Rather 
than growing in amplitude as dictated by the lowest-order 
energy, the wrinkled state collapses into a fold. 

The enhanced localization arises because of the unex- 
pected importance of non-leading energetic terms at in- 
definitely small A. An extended state, however small its 
amplitude, can be rearranged to produce a localized state 
of nonvanishing amplitude, for which nonleading contri- 
butions to the energy are significant. (As we noted, any 
uniform wrinkle must have vanishing amplitude and hence 
vanishing anharmonic energy.) 

The localized buckling studied here resembles the buck- 
ling of a thick elastic film that is put under compressive 
strain 20,21 . In this case as in ours, linear response sug- 
gests an undulating surface. However, higher-order parts 
of the energy lead to a localized crease of nonzero slope but 
arbitrarily small depth. An analogous localized instability 
occurs in an elastic solid under tensile stress. Typical mod- 
els of elastic solids are linearly unstable to the formation of 
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cavities of arbitrarily small size [22 . Our case differs from 
these two in that our localized state occupies a region of 
diverging size as the threshold is approached. This "pro- 
gressive localization" bears closer resemblance to that seen 
in the nonlinear Schrodinger equation for wave- function ip: 

Here any negative energy E produces a localized bound 
state with a spatial extent that diverges as E — > (albeit 
with no oscillation) |23] , 



Conclusion 

The incipient folding instability described above opens 
several promising directions. The energy and incipient shape 
found here form the basis for predicting the time depen- 
dence of folding. The observed abrupt appearance and ar- 
rest of folds in time has long been a puzzle [5] . The unex- 
pected simplicity of the energy versus displacement sug- 
gests an underlying symmetry of the system as yet undis- 
covered. Our work on these problems is in progress. 

The authors warmly acknowledge insightful discussions with 
Ka Yee Lee, Luka Pocivavsek, Jin Wang, Enrique Cerda, and 
Benny Davidovitch. This work was completed at the Aspen 
Center for Physics. It was supported in part by the US-Israel 
Binational Science Foundation under Grant Number 2006076, 
and in part by the National Science Foundation's MRSEC Pro- 
gram under Award Number DMR 0820054. 



infinite one; (b) we expand the expressions in small k but 
keep terms that include kL intact. The resulting energy is 

E = 2(l + 2 K 2 )A-^KCoth(^^j A 2 . (A.2) 

Minimization of E with respect to k yields the follow- 
ing equation: 

AL „. r . . sinhx — x , . „. 

nL=—f{KL), /(«) = _ E _ I . (A.3) 

For A < A c — 48 /L this equation has only the trivial 
solution, k = 0, corresponding to the wrinkled state. For 
A > A c a solution of finite k appears, increasing contin- 
uously from zero and corresponding to the folded state. 
Thus, sufficiently close to the transition we may expand 
eq. (|A.2[) in small k to get 

E ~ E + 4Z\(1 - A/A c )n 2 + {L 3 A 2 /720)K i , (A.4) 

where Eq = 2A — A 2 /L. This Landau-like energy makes 
evident the second-order nature of the transition, k serv- 
ing as the order parameter. From eq. (|A.4[) we obtain 

JO, A< A c 

K ~ \ L^A-^iA - 4c) 1 / 2 , A>A C 

E ~ E °~ \-L- l (A-A c ) 2 , A>aI (A ' 5) 

At A = A c there is a discontinuous jump in d 2 E/dA 2 = 
dP/dA 



Appendix A. Finite-size effect 

The analyses in sects. [3] and 0] have been restricted to the 
limit of a very large system size, 1 <C n^ 1 <C L. However, 
above and sufficiently close to the wrinkle-to-fold transi- 
tion (i.e., for A > A c = 48/L), n is arbitrarily small, and 
the assumption kL 3> 1 must break down. We now treat 
this limit of very small displacement, where (k _1 ,_L) 3> 1 
but nL may be finite. 

For such minute deformations an 'evanescent wave' 
Ansatz, 

h(s) = Asin(fcs)e- K|s| , (A.l) 

gives the correct qualitative behavior as regards the scal- 
ing laws while simplifying the required integrations. An 
odd function has been chosen in eq. (|A.1|) to avoid a diver- 
gent curvature at s = 0. We set k = k c = 1 and L = 27m, 
where n is a large integer (i.e., the system size is compati- 
ble with a large integer number of wrinkling wavelengths). 
We repeat the same procedure from sect.H] — substitute 
eq. (|A.1[) in eq. @, calculate E and A to second order in 
A 2 , express A 2 as a function of A, and substitute it back 
in the energy to get E(A). There are two differences in the 
current calculation: (a) the integrations in eq. ^ are per- 
formed over the finite interval (— L/2, L/2) rather than an 
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